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Abstract 

The recent experimental data obtained by the OBELIX group on pD and 
p'^He total annihilation cross sections are analyzed. The combined analysis of 
these data with existing antiprotonic atom data allows, for the first time, the 
imaginary parts of the S-wave scattering lengths for the two nuclei to be ex- 
tracted. The obtained values are: Im a^'^ = [— 0.62±0.02(stat) ±0.04(sys)] fm 
for pD and Im Oq'^ = [—0.36 ± 0.03(stat)to;n(sys)] fm for p^He. This analysis 
indicates an unexpected behaviour of the imaginary part of the p-nucleus S- 
wave scattering length as a function of the atomic weight A: 



Im Oq^I (pp) > |Im Uq'^I (pD) > |Im Og'^l (p'^He). 



PACS numbers: 25.43.+t, 36.10-k 
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1 Introduction 



Recent data obtained in the experiments on pD and p^He annihilation in flight at low 
energy and on the measurement of the shift and the width of the 15* level of the pD 
atom gave a first indication of a quite unusual phenomenon. Indeed, the imaginary 
part of the pD scattering length appears to be smaller than the pp one . The performed 
experiments were the first investigations of the pD system at low energy; their realization 
was a very difficult task and the precision and accuracy obtained were limited. 

The first aim of this work is to combine the data coming from these experiments 
in order to extract improved information on the imaginary part of the pD scattering 
length. As we have shown in our previous paper [Q, devoted to the pp annihilation 
at low energies 0, H, a successful analysis can be done within the scattering length 
approximation written for the systems with Coulomb attraction. These two approaches, 
in-fiight annihilation and antiprotonic atom experiments, were shown to give coherent 
results and can be used as complementary. The second aim of this work is to perform a 
first phenomenological analysis of the antiproton annihilation on heavier nuclei in S-wave 
and in P-wave. 

The values of the pD and p^He annihilation cross sections used in this paper are 
reported in Tab. These are the only p-nucleus annihilation data available at incident 
momentum below 200 MeV/c. The data concerning shifts and widths of antiprotonic 
atoms will be reported in the text. 

Table 1: Values of the pD and p^He total annihilation cross sections used in this work, 
multiplied by the square of the incoming beam velocity (3 for different p incident momenta. 
In the data from [|l|, in addition to the quoted statistical and systematic errors, an overall 
normalization error of 2.5% has to be considered. 



gaseous 
target 
(ref.) 


p incident 
momentum 
(MeV/c) 


ann 

(mbarn) 


D2 i 


69.6±1.5 
45.7±3.5 
36.3±5.1 


3.45±0.08 (stat) ±0.15 (sys) 
2.12±0.06 (stat) ±0.33 (sys) 
1.96±0.08 (stat) ±0.55 (sys) 


'Hei 


70.4±1.3 
47.0±3.3 


4.63±0.10 (stat)±0.19 (sys) 
2.45±0.10 (stat)±0.35 (sys) 


'He 


45.0±5.0 


3.1±0.7 



The structure of the article is the following. In section 2 we present all the necessary 
scattering length formalism for systems with Coulomb attraction. In section 3 the last 
experimental data on the pp annihilation cross sections are analyzed and the pp low 
energy parameters, extracted both from in-flight annihilation and from atomic data, are 
shown to be in excellent agreement. Section 4 is devoted to the combined analysis of 
the pD data coming from in flight annihilation and atomic experiments; in this section 
the imaginary part of the pD scattering length is obtained. In section 5 we apply the 
analogous procedure to the p^He annihilation data and we obtain the imaginary part of 
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the p^He scattering length. In addition, we examine the behaviour of p-nucleus scattering 
parameters as a function of the atomic weight A, for S-wave and P-wave, for different 
nuclei. Finally, the conclusions contain a brief summary of the results. 



2 Scattering length approximation 

The starting point to develop the scattering length approximation is the relation between 
the X-matrix for a given orbital momentum / and the strong interaction phase shift in 
presence of Coulomb forces 0: 

j^^^ = gi{v)q''^'Kiv)cot6r - 2vh{v)], (1) 

where q is the center-of-mass momentum, B the Bohr radius, 77 = 1/ qB 

9o{v) = 1, 

giiv) = n + ^ = 1'2' 

mil V ^ J 



with the digamma function The i^-matrix is related to the S'-matrix by 

^ l + igi{r])q^^+^w{r]yK 
l-igiiri)q'''+^w{ri)K 

with w{ri) = Cq^t]) + 2irjh{vi). 

The scattering length approximation used at low energies is equivalent to the replace- 
ment of the ii'-matrix by a constant: 

Within this approximation, the annihilation cross section a^„„ for a given orbital momen- 



tum / takes the form |TU[: 



q ^ann = i^l + 1)4. 1 1 _ (^),.e| 2 ' (2) 

The scattering length approximation allows the low energy parameters to be obtained 
also from atomic data. Following Trueman [|ll|], it is necessary to replace in (|^) q by 
—i\/2iiE, fi being the reduced mass, E = Eni + AEni the exact position of the Coulomb 
level Eni shifted and broadened by the strong interaction AEni, and cot5f^ by i: 

-\-, = -9liv)q''^'[C'oiv)^~2vh{v)]. (3) 

Within this approximation (when one neglects effective range corrections), this expression 
is exact. If one supposes that AEni/E^i <^ 1 or, equivalently, a^^/i?^'+^ ^ 1 one obtains 
different approximate relations usually called Deser [0 (for S-wave in the first order) or 
Trueman [TT[] (for S- and P-wave in higher orders) formulas. 
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3 pp annihilation 

To illustrate the agreement between the values of the low energy parameters extracted 
from the in-flight annihilation and from the atomic data, let's start from the pp system, 
where the experimental data at very low energy are more abundant and precise. This 
procedure was already applied to the pp system and is described in detail in 0]. In the 
present analysis, we add recently measured experimental points ||13[. The results of the 
fit are presented in Fig. |l|. As in the previous analysis [Q, the experimental point at 43.6 
MeV/c, which suffered from possible unknown systematics |T^, was not used for the fit. 




Figure 1: Values of the total pp annihilation cross section multiplied by the square of the 

(ir). The error 
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incoming beam velocity. Experimental data are from [g] (a), |§] (■ 

bars represent the quadratic addition of the statistical error and of the systematic error 
interval divided by \/T2. Moreover the data are affected by an overall normalization error: 
3.4% for the data from [^, ^] and 2.5% for the data from |jl3|. The theoretical curves are 



the result of the present work: the full line is the total annihilation cross section, 
dashed line represents the S-wave contribution. 



the 



The fit was performed by means of the MINUIT program |14 and provided the fol 



lowing best fit values for the imaginary parts of the pp scattering length (S-wave) and of 
the scattering volume (P-wave): 

Im a'o" = -[0.69 ± O.Ol(stat) ± 0.03(sys)] fm; 

(4) 

Im af = -[0.75 ± 0.05(stat) ± 0.04(sys)] fm^ 

with the value of = 0.25 per point. In the fitting procedure only the statistical errors 
were accounted for; their propagation produced the errors quoted as statistical in the 
values of the best fit parameters. The errors quoted as systematic come from the overall 
normalization error of the experimental data. Let us remind the reader that the following 
additional hypotheses were used in the fitting procedure: the parameters are the spin 
averaged ones, | Re Og'^l = | Im a^'^l and the real part of the scattering volume can be 
neglected. 

It is very instructive to compare the values obtained for the two parameters with the 
ones obtained from the atomic data. The last world averaged value for the shift and the 



width of the IS* atomic level of the pp atom is [16 



AEis + = [-(0.721 ± 0.014) + i(0.548 ± 0.021)] keV 

which gives, by means of the Trueman formula, the imaginary part of the scattering 
length: 

Im a'o" = -(0.694 ± 0.027) fm, 

in excellent agreement with the result obtained from the in-flight annihilation data (^. 

The imaginary part of the P-wave scattering volume can be extracted from the width 
of the 2P atomic level of the pp atom. The majority of the atomic experiments obtain 
this value indirectly through the intensity balance procedure which gives, actually, the 
lower limit for this parameter (see discussion in |TB|, |T^). The world averaged value 
obtained by this method is fir 



T2P = (32.5 ±2.1) meV 

which corresponds to: 

Im af = -(0.66 ± 0.04) fm^ 

This value is smaller than the value obtained from the in-fiight annihilation data (^). 

However, the unique direct measurement of the width of the 2P atomic level, which 
has been recently performed |T^, gives: 



Fap = (38.0 ± 2.8) meV 

which corresponds to: 

Im af = -(0.77 ± 0.06) fm^ 

in excellent agreement with (^. 

The main conclusion of this analysis is that the results obtained from these different 
experimental approaches are in quite good agreement; therefore these experiments can be 
considered as complementary. 
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4 pD annihilation 



The first measurement of tlie IS* level of tlie pD atom p| gave a very surprising result: 
the width of this level seems to have the same size as the corresponding one of the pp 
atom. The shift and the width of the IS atomic level were found to be: 

AEis + = [-(1.05 ± 0.25) + z(0.55 ± 0.37)] keV, 

which give the following scattering length for the pD system: 

ag^ = [(0.7 ± 0.2) - i{OA ± 0.3)] fm. 

Unfortunately, the experimental errors are too large to allow for an unambiguous com- 
parison with the pp scattering length. 

The precision of the imaginary part of the scattering length can be improved signifi- 
cantly if the information coming from atomic measurements are combined with the data 
on pD annihilation in flight. The general idea is quite simple: the P-wave parameters 
are fixed from the atomic data (which are quite precise in this case), the real part of the 
scattering length is taken from the atomic data too (the annihilation cross section is not 
very sensitive to this parameter). Thus we can perform a fit to the annihilation cross 
section data with only one free parameter: the imaginary part of the S-wave scattering 
length. 

The shift and the width of the 2P atomic level were measured directly with approxi- 
mately ten percent accuracy [ P^ : 

AE2P + = [-(243 ± 26) + i(245 ± 15)] meV 

which correspond to the following scattering volume: 

af = [(3.3 ± 0.3) - i(3.18 ± 0.18)] fml 

The results of the fitting procedure are presented in Fig. the procedure provided 
the following best fit value for the imaginary part of the scattering length: 

Im ag^ = -[0.62 ± 0.02(stat) ± 0.05(sys)] fm 

with a rather large value of the x^? which amounts to 8.3 per point. 

In the fitting procedure only the statistical errors were accounted for; their propagation 
produced the error quoted as statistical in the value of the best fit parameter. The quoted 
systematic error originates from two main sources of systematics. The first one is the error 
on the determination of the imaginary part of the scattering volume (giving 0.03 fm) and 
the second one is connected to the overall normalization error of the data (giving 0.04 
fm). These two errors are added quadratically. On the contrary, the imaginary part of 
the scattering length is practically insensitive to the variation of its real part, within the 
experimental errors taken from the atomic experiment. 

The agreement of the theoretical curve with the experimental data, within the statisti- 
cal errors, is not very good, as shown by the large value of the x^- This discrepancy could 
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Figure 2: Values of the total pD annihilation cross section multiplied by the square of the 
incoming beam velocity. Experimental data are from (•). The error bars represent the 
quadratic addition of the statistical error and of the systematic error multiplied by \fvi. 
Moreover the data are affected by an overall normalization error of 2.5%. The theoretical 
curves are the result of the present work: the full line is the total annihilation cross 
section, the dashed line represents the S-wave contribution, the dotted line the P-wave 
contribution. 



be due to at least two reasons. The experimental points, especially at the lowest values 
of the incident momentum, are affected by large individual systematic errors, due to the 
difficulty of separating annihilations coming from the different momentum components of 
the antiproton beam Moreover, at the lowest value of the incident momentum, the 
spread of the projectile momentum distribution is quite large (see Tab. P; as the annihi- 
lation cross section, in this momentum range, is rapidly increasing with the lowering of 
the incident momentum, the measured value of the cross section should be considered as 
an average value over the momentum interval, rather than the cross section value at the 
center of the interval. 

In the error bars reported in Fig. § the two effects previously mentioned are acconted 
for and the theoretical predictions appear in fair agreement with the experimental data. 
Therefore, in spite of the large value of the of the fit, we consider that the result 
obtained for the best fit parameter could be confidently accepted. 
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In conclusion, the combined analysis of the in-flight annihilation data and atomic data 
allows the knowledge of the imaginary part of the pD scattering length to be improved 
significantly. Moreover, the indication obtained from the atomic experiments that this 
value does not exceed the pp one is confirmed. 

As a final comment, let us remark that, within a naive geometrical approach to anni- 
hilation, one could expect the imaginary part of the pD scattering length to be approxi- 
mately equal to the sum of the pn and pp scattering lengths: 

Im ao(pD) ^ Im ao(pn) + Im ao(pp). 

These results show that this vision is completely wrong especially if one takes into account 
the quite large value of the imaginary part of the rip (or equivalently pn) scattering length 

Im ao(np) = -[0.83 ± 0.07(stat)] fm, 

which was obtained from the data on lip annihilation |TH . 



5 p He annihilation and P-wave in different nuclei 

The ground state of the p^He atom is experimentally unaccessible. Therefore one cannot 
obtain any information about the S-wave scattering parameters from atomic experiments. 
Nevertheless, it is possible to evaluate the imaginary part of the p^'He scattering length 
if one combines the atomic information with the data on in-flight annihilation []I[ . 
The P-wave scattering volume: 

af = [-(3.4 ± 0.4) - i(4.4 ± 0.5)] fm^ 

as well as the imaginary part of the D-wave scattering parameter : 

Im a^" = -(1.42 ± 0.06) fm^ 

can be extracted from the atomic data [jl9[. Unfortunately, there is no experimental 
information about p^He interaction in S-wave. To perform the fit, we choose the value 
of the real part of the S-wave scattering length Re Og'^ = 1.0 fm, which is only slightly 
higher than the corresponding parameter for the pp and pD systems. Given the large 
arbitrariness on this parameter, we assume an error of 50% on its value (±0.5 fm). 

Thus, we can perform a fit of the available p^He annihilation cross section data with 
only one free parameter: the imaginary part of the scattering length. The result of the 
fit is presented in Fig. the experimental data considered are from [^ and from [|^, a 
previous low statistical measurement of the p'^He annihilation cross section performed at 
LEAR with a streamer chamber. 

The fit to the annihilation cross section data provided the following best fit value for 
the imaginary part of the scattering length: 

Im a^^ = [-0.36 ± 0.03(stat)+°:i?(sys)] fm. 

with a value of x^=l-8 per point. In the fitting procedure only the statistical errors 
were accounted for; as in the previous cases their propagation produced the error quoted 
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Figure 3: Values of the total p^He annihilation cross section multiplied by the square 
of the incoming beam velocity. Experimental data are from [|l| (■) and (a). The 
theoretical curves are the result of the present work: the full line is the total annihilation 
cross section, the dashed line represents the S-wave contribution, the dotted line the 
P-wave contribution, the dashed-dotted line the D-wave contribution. 



as statistical in the value of the fit parameter. Here the systematic error contains two 
contributions. The first comes from the error on the imaginary part of the scattering 
volume and from the uncertainty on the real part of the scattering length; the second 
comes from the overall normalization error of the data. The contribution coming from 
the imaginary part of the D-wave scattering parameter is negligible. As a final remark, 
let us emphasize that this is the first experimental evaluation of the imaginary part of the 
p'^He S-wave scattering length. 

Figure ^ summarizes the values of the imaginary part of the p-nucleus scattering 
length as a function of the atomic weight A. This function seems to be a decreasing one. 

It is interesting to compare this unusual behaviour with the behaviour of the imaginary 
part of the p-nucleus P-wave scattering volume as a function of the atomic weight. For 
pp, the information is obtained from atomic and annihilation experiments, as described 
in section 3. For heavier nuclei (D, He, Li) the data come from the atomic experiments 
16, 19, pOi, where the shifts and the widths of the 2P levels are measured directly. The 



experimental widths as well as the imaginary parts of the scattering volumes, calculated 
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Figure 4: Absolute values of the imaginary part of the antiproton-nucleus scattering 
length as a function of atomic weight A. Statistical and systematic errors are added 
quadratically. 



through the first order of the Trueman formula, are presented in Tab. 



Table 2: p-nucleus 2P level widths and Im al'^ calculated in the first order of the Trueman 
formula. 



System 


ref. 


^2p 

(meV) 


— Im a^^ 
(fm^) 


pH 
pD 
p^He 
p^He 
p^Li 
p^Li 


[ 
[ 
[ 


16 
16 

19 
20 
20 




38.0 ±2.8 
489 ± 30 
(25 ±9) ■ 10=^ 
(45 ±5) ■ 10=^ 
(444 ±210) ■ 10^ 
(456 ± 190) ■ 10^ 


0.77 ±0.06 
3.18 ±0.20 
3.1 ± 1.1 
4.4 ±0.5 
4.3 ±2.0 
4.1 ± 1.7 



These values, as a function of the atomic weight, are presented in Fig. ^. Unfortunately, 
the large experimental errors do not allow for an unambiguous conclusion about the 
form of the functional dependence. Nevertheless, this function seems to have, at least, a 
saturation-like behaviour. 

An analogous question can be put forward for the D-wave. The experimental atomic 



information for heavy nuclei is quite abundant (see for review ||2T|). However, the scatter- 
ing length approximation used here to extract the scattering parameter af^ is no longer 
valid. When one works quite far from the threshold, it is necessary to take into account 
higher terms in the development of the i^'-matrix and to use an effective range approxi- 
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Figure 5: Absolute values of the imaginary part of the antiproton-nucleus scattering 
volume as a function of the atomic weight A. 



mation instead of the scattering length one. Thus, the K-matrix must be written as: 



Here q is the momentum of the corresponding Coulomb level: q — iqs — i/nB. For light 
nuclei, like hydrogen or helium, for which B = 57.6 fm and B = 18 fm respectively, a 
correction coming from the second term is negligible, if one supposes that rf^ is of the 
order of 1.0 fm. For heavier nuclei, for instance ^'^F, B = 3.37 fm and this value is of the 
same order as all other parameters in this expression. Therefore, the development has 
no sense and the scattering parameters cannot be extracted from the results of atomic 
experiments. Note that the experiments measuring annihilation cross sections have no 
such problem because, in principle, any value of q can be chosen. 



The proposed combined analysis of the recent p-light nucleus annihilation data and an- 
tiprotonic atom data allows our knowledge of the low energy parameters in these systems 
to be improved. The imaginary part of the spin-averaged pD scattering length can be 
obtained: 




6 Conclusions 



Im a^"(pD) 



[0.62 ± 0.02(stat) ± 0.05(sys)] fm. 
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Moreover, for the first time, it becomes possible to evaluate the imaginary part of the 
spin-averaged p^He scattering length: 

Im ag"(p^He) = [-0.36 ± 0.03(stat)lgi?(sys)] fm. 

Compared to the the imaginary part of the spin- averaged pp scattering length 

Im a^^(pp) = -[0.69 ± O.Ol(stat) ± 0.03(sys)] fm, 

obtained from annihilation data and in agreement with antiprotonic atom experiments, 
these values indicate the presence of a quite unexpected phenomenon: the absolute value 
of the scattering length seems to be a decreasing function of the atomic weight. 

A naive geometrical vision of p-nucleus annihilation would suggest a value of the p- 
nucleus scattering length increasing with the nucleus size. This naive picture is also wrong 
when one analyses the imaginary part of the P-wave scattering volume: the function seems 
to have, at least, a saturation-like behaviour. 

To confirm this unusual phenomenon involving low energy p-nucleus interactions, it 
would be necessary to perform new measurements with higher statistics. This experimen- 
tal information, expanded to heavier nuclei, could be obtained with the new AD facility 
at CERN. 
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